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In this paper we calculate the block entanglement entropies of spin models whose ground states 
have perfect antiferromagnetic or ferromagnetic long-range order. In the latter case the definition 
of entanglement entropy is extended to properly take into account the ground state degeneracy. We 
find in both cases the entropy grows logarithmically with the block size. Implication of our results 
on states with general long-range order will be discussed. 



I. INTRODUCTION 

Entanglement is the hallmark as well as most counter- 
intuitive feature of quantum mechanics. Among various 
ways to quantify entanglement, bipartite block entangle- 
ment entropy emerged as a concept of central importance 
in quantum information science [ij, and has been receiv- 
ing growing attention in other branches of physics re- 
cently. For example, it was suggested to be a possible 
source of black hole entropy @> La]. In condensed mat- 
ter or many-body physics, the entanglement entropy has 
been increasingly used as a very useful and in some cases 
indispensable way to characterize phases and phase tran- 
sitions, especially in strongly correlated phases [3]. In this 
context the most important result is perhaps the so-called 
area lawHHl], which states that in the thermodynamic 
limit, the entropy should be proportional to the area of 
the boundary that divides the system into two blocks. 
There are a few very important examples [1,0,0] in which 
the area law is violated, most of which involve quantum 
criticalityQ; the specific manner with which the area law 
is violated is tied to certain universal properties of the 
phase or critical point. In some other cases, important 
information about the phase can be revealed by studying 
the leading correction to the area law@,H,[§]; for example 
this is the case for topologically ordered phases @, [!]• 

Comparatively there have been relatively few studies 
of the behavior of entanglement entropy in states with 
traditional long-range order [13 US [l2| . This is perhaps 
because of the expectation that ordered states can be well 
described by mean-field theory, and in mean-field theory 
the states reduce to simple product states that have no 
entanglement. In particular in the limit of perfect long- 
range order the mean-field theory becomes "exact" , and 
the entanglement entropy should vanish. In this paper 
we will show that this is not the case, and interesting 
entanglement exists in states with perfect long-range or- 
der. We will study two exactly solvable spin-1/2 models: 

(i) An unfrustratcd antifcrromagnct with infinite range 
(or constant) antiferromagnetic (AFM) interaction be- 
tween spins in opposite sublattices, and ferromagnetic 
(FM) interaction between spins in the same sublattice; 

(ii) An ordinary spin-1/2 ferromagnet with arbitrary FM 
interaction among the spins [Hj]. While the ground states 
have perfect long-range order for both models, we show 
that they both have non-zero entanglement entropy that 



grow logarithmically with the size of the subsystem. 

The paper is organized as follows. In Sec. II we intro- 
duce the model (i) and present its exact solution. In Sec. 

III we calculate the reduced density matrix of a subsys- 
tem to obtain the exact expression of the entanglement 
entropy for model (i). Based on those results, in Sec. 

IV we study the scaling behavior of the entropy as the 
system size tends to infinity under two different partition 
limits and compare it with numerical calculations. Sec- 
tion V is devoted to model (ii). Finally, in Sec. VI we 
summarize and discuss the results of this paper. Some 
mathematical definitions, notations and details are given 
in the appendices. 



II. ANTIFERROMAGNETIC SPIN MODEL 
AND GROUND STATE 

We consider a lattice model composed of two sublat- 
tices interpenetrating each other as in Fig. [TJ with inter- 
action of infinite range, i.e., every spin interacts with all 
the other spins in the system, with interaction strength 
independent of the distance between the spins. Within 
each sublattice, the interaction is ferromagnetic, and be- 
tween the sublattices the interaction is antiferromagnetic; 
as a result there is no frustration. The Hamiltonian is 
written as, 

H = —Ja ^2 Si-Sj — Js ^ Si-Sj + Jo ^2 

i.jeA i,j&B i£A.j£B 

(1) 

with Ja,Jb,Jo > 0. The ground state of Eq. |T]) can be 
solved in the following manner [l^]. Let 



Sa = ^2 S tl 

i£A 



Sb — ~^2^i, S — Sa + Sb, 



ieB 



then the Hamiltonian can be written as, 
H = —JaS a — JbS b + JqSa ■ Sb 

= - J A S 2 A - JbS 2 b + -^{S 2 - s a - s b)- 

Since [S,H] = [S A ,H] = [S B ,H] = 0, eigenstates of H 
can be chosen to be simultaneous eigenstates of S 2 , S A 
and S B , with angular momentum eigenvalues S, Sa and 
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FIG. 1: Two-sublattice model: two sublattices labeled A and 
B interpenetrating each other. 



Sb respectively, such that the energy eigenvalue is, 



^S(S+1)-(J A +^)S A (S A - 



1)-(Jb + y)Sb(S b + 1). 



(3) 

To minimize the energy, we must first maximize S A and 
Sb, and then minimize S, i.e., S = \S A — Sb\,S a = 
^,Sb = -sr. Here N A and Nb represent the number 
of corresponding sublattice sites. Physically this means 
that spins in the same sublattice are all parallel to each 
other, while spins in opposite sublattices are antiparallel. 
Thus the ground state is, 



\S A S B ; Sm) 
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(4) 



In this paper we only consider the simplest case with 
N A = Nb = N, thus the total system size is 2N. 
Then the ground state is reduced to an antiferromagnetic 
ground state which has zero total spin, 



\S a Sb\ Sm) 



N N 



;00). 



This ground state has perfect Neel order, as manifested 
by the spin-spin correlation function, 



(Si ■ Sj) = -, i,j 6 A or i,j 6 B: 

(Si ■ Sj) = _i _ -L - -1 i 6 A and j e B. 



(6) 



III. THE REDUCED DENSITY MATRIX AND 
ENTANGLEMENT ENTROPY 

We divide the system spatially into two subsystems 
which are labeled 1 and 2 respectively, as shown in Fig. 
[21 and study the ground state entanglement entropy E 
between these two subsystems. E is defined to be the 
von Neumann entropy of either of the subsystems, E\ 
or E2, which can be calculated from the reduced density 
matrix, 



E = Ei — E2 



-tr(pi lnpi) = -tr(p 2 lnp 2 ), (7) 



where p\ is the reduced density matrix of the sub- 
system 1, obtained from the density matrix p = 
IyT' 00)(tT>00| °f the whole system by tracing out 
degrees of freedom of the other subsystem, p\ = tr^ 2 
and vice versa. 



ABAB A B 

Subsystem 1 



t I ~t I t I 

A B ABAB 

Subsystem 2 



FIG. 2: Bipartite division of the system. We divide the system 
spatially in to two parts and label them subsystem 1 and 
subsystem 2, respectively. One of the main tasks of this paper 
is to evaluate the entanglement entropy between these two 
parts. 



To solve for the explicit form of the reduced density 
matrix, we proceed as follows. First, we further decom- 
pose the system into four parts, S Al , S A2 , Sb ± , Sb 2 , 
with 



= DieBAJSl &B 2 = J2ieBAi<£2 



a.. 



(8) 



Therefore, these operators satisfy the following relations, 



S A = S Al + Sa 2 
Si = S Al + Sb! 



Sb = Sb! 
S2 = Sa 2 



f Sb 2 
Sb 2 - 



(9) 



(5) Here we note that, as discussed in the previous section, 
the spin state within each sublattice is ferromagnetic. 
This means that not only must the total spin quantum 
numbers of S A and Sb take their maximum values, but 
the total spin quantum numbers of S Al , S A2 , Sbx and 
Sb 2 must also take their maximum values. More im- 
portantly, these values are thus fixed, which enables us 
to treat the operators S Al , S A - 2 , Sbx and Sb 2 as four 
single spins, and in what follows we shall denote these 
operators by their corresponding spin quantum numbers. 
The problem is then that we are given a four spin state 
in which the spins S Al and S A2 are combined into a state 
with total spin S A and the spins Sb ± and Sb 2 are com- 
bined in a state with total spin Sb, and then these two 
states are combined into a total singlet (resulting in the 
ground state of our long-range AFM model), and we must 
express this state in a basis in which the spins S Al and 
Sb ± are combined into a state with total spin Si and S A2 
and Sb 2 are combined into a state with total spin S^. 
This change of basis involves the familiar LS-jj coupling 
scheme. 

We proceed by first rewriting the density matrix de- 
scribing the ground state in terms of the basis states 
ISimiiS^m^) by appropriate insertions of the unit opera- 
tor 

(p), \Sim 1 S 2 m 2 )(S 1 
S17711S2T12 



N N N N 
p=\S A S B ;Sm)(S A S B ;Sm\ = | — — ; 00)< — — ; 00| 

= ^ \Sim 1 S2m2)(SimiS2m 2 \SASB; Sm) 

Simi 
S 2 m 2 

x (S A S B ;Sm\ £ \S'im'iS' 2 m' 2 ) {SW^ 

S' 2 m 2 

^2 X s imi s 2m2 ys[rn[s' 2 7n'JSim 1 S2m2)(S' 1 m' 1 S! 2 rn' 2 

(10) 



A Simi =(-l) Sl - m V(25i + 1)(2S A + 1)(2S B + 1) 



Sa! S Bl Si 
x | Sa 2 S b , 2 Si 
Sa S b S 



(14) 



Sim 1 S 2 m 2 
S[m' 1 S 2 m 2 



where As lTOl s 2TO2 = {SimiS 2 m 2 \SAS B ; Sm) . We then 
insert another unit operator 

\S 1 S 2 ;S Q M)(S 1 S 2 ;S M\ 

to express this matrix element as the product of 
a Clebsch-Gordan coefficient and an LS-jj coupling 
coefficient [l6|. 

{SimiS 2 m 2 \SAS B ] S m) 

= (SiniiS2m2\SiS2SoM)(SiS2SoM\SAS B Sm) 



, In our case S = 0, so the 9-j symbol can be expressed 
more simply in terms of a 6-j svmbolfl8j (see also Ap- 
pendix B), 

'a b c\ 

d e / =5 cf S gh (-iy +d+c+9 ((2c+l)r? 
\9 h 0/ ( 15 ) 

a b c 



x (20 + 1)" 



e d .9 



So,M 



J So S mi +7712)771^ M,m 



(SimiS2m2\SiS2Sm) 



Sai S Bl Si 
Sa 2 S B2 S 2 
Sa S b S 

(11) 



In our particular case, within each sublattice, the state 
is ferromagnetic, Sa = ^j, S B = ^j. Let 2N\ be the size 
(number of lattice sites) of subsystem 1. Without losing 
any generality, we can let N\ ^ N. Applying Eq. JTJ 
we then have, 

/ jV-JVi jV-TVi c 

\ Sl , mi = (-i) Sl - mi VW+i~)( A A I 

\ 2 2 2 

(1§ 

Using a symmetry property of the Racah coefficients fl7l| 
W(abcd; ef) = W(acbd; fe), we have, 



' N-N x _/Vi jV 

\ Simi = (-i) Sl - mi V(NTT){ A 2 

2 2 



SUj &i Si 

Here SU 2 Sb 2 S2 
Sa S b S 
X-coefficient defined as 



is the LS-jj coupling coefficient or 



Then we employ the following relation [171]: 

(2a)!(26)! 



(17) 



W(abcd, a + b, /) 



h si ji 
h s 2 j 2 
L S J 



(lll 2 SlS2LSJM\liSll2S2jlj2JM). (12) 



(2a + 2b + l)!(c + d - a - b)\ 
(a + b + c- d)l(a + b + d~ c)!! 
(a + f-c)\(a + c + f + l)\(b + d-f)\ 

(c + f-a)\(d + f-b){a + b + c + d+l)l 



Since S = m = 0, we have Si = S 2 , mi 



-m 2 . So 



(b + f-d)\(b + f + d+l)\(a + c-f)\ 



{Sinn; S2m 2 \SiS2] Sm) 



(-iy- 



and for simplicity, 



(18) 



V252+1 

we can now suppress the S2 and m2 indices and repre- 
sent As imi 5 2m2 by Xsmn without causing any ambiguity. 
Therefore, 



This gives 
A Sl , mi = (-l) s — 



(N + 1)(N - NiJlN^. 



*Simi 



(-1) 



Si— mi 



V25T+T 



S Al S Bl Si 
Sa 2 S B2 Si 
Sa S b S 



(13) 



(N-Ni- Si)\(N -N1 + S1 + 1)! 
(N — Ni)\Ni\ 



(Ni-Si)\(Ni + S 1 + l)\ 



(19) 



In the following calculation, for consistency, we adopt the 
convention of Wigner 6-j and 9-j symbols, the definitions 
of which, and the relation to the Racah coefficients and 
LS-jj coupling coefficients (or X-coefficients) , are given in 
Appendix A. 

Following Wigner's convention, Eq. JTSj) can be written 



Now with the expression above, we can directly trace 
out the degrees of freedom in subsystem 2 from Eq. (fTU|) . 
with the result 



Pi 



= tr {2 ){p)= ^2 |As 1 m 1 | 2 |5'imi)(S'iTOi| 



(20) 



Simi 



4 



The bipartite entanglement entropy between subsystems 
1 and 2 is then given by 

E = Ek = -Y, I**«h | 2 ln(|A Simi | 2 ). (21) 

Si ,mi 

Here we note that, although As imi is written with an 
explicit mi dependence, the actual expression is indepen- 
dent of mi. As a result, we can eliminate the summation 
over mi from (|2ip by multiplying by a factor of 25*1 + 1. 
The final expression for the entanglement entropy is then 

E = E\ = — ^(25i + l)|A Simi | 2 ln(|A Simi | 2 ). (22) 

Si 

In the following section we will first derive the asymp- 
totic behavior of the entanglement entropy E in certain 
limits using the exact results above, and then present 
the results of our numerical calculation of E with finite 
system sizes. 



IV. ASYMPTOTIC BEHAVIOR AND 
ENTANGLEMENT ENTROPY 



this criteria is satisfied quite well in both cases. Then by 
the normalization condition 



Y^A{2Si + l)e 



-sli 



2 r 1 Li' 

N -N! ^ N± ■ 



Si 



A(2Si + l)e~ Sl{w ^ +1 ^ ) dS 1 = 1 



1 / w(N-Ni)Ni • 



\ 2 

A Simi 



sf sf 



The entanglement entropy is therefore given by 

Ei~- [ A(2S 1 + i) e - s i ( ^r + ^r ) 
Jo 

x In ( J 4e~ 5?( ^r + ^ ) ) dSi 



(24) 



we obtain the (inverse) prccoefficient jjj- ~ (Ni — ^> 
-jj ). Therefore, 



-e ™i "-"i . (25) 



(26) 



In this section we consider two limiting cases: 

(i) A^ = N<i = N, this case gives the saturated entropy 
at fixed N since intuitively E should increase with the 
subsystem size. 

(ii) 1 <C iVi <C N, in this limit we are considering sys- 
tem's entanglement with its (much larger) environment, 
and gencrically we should be able to find that the en- 
tropy should be independent of the total system size as 
TV — > oo, which is indeed what we find. 

To obtain the asymptotic behavior of the entangle- 
ment entropy, we first consider the asymptotic behavior 
of As imi , then in the large N limit. Then the summation 
over Si which runs from to N/2 can be replaced by an 
integral. As we will see, the distribution of Xs imi is pro- 
portional to a Gaussian function with respect to Si , thus 
the bounds of integration can be extended to from to 
+oo. 

Our concern is the distribution of As imi with respect 
to Si, therefore, we extract the dependence on Si and 
then use the normalization condition y% m Ac m = 
J2s + l)A| imi = 1 to obtain the normalization con- 
stant. From Eq. (|T9|) we can write, 



\ 2 

A Simi 



(2iVi)! 



(JVi-5i)!(JV + Si + l)! 

(2N-2N{)\ 
(N-Ni- Si)l(N -N1 + S1 + 1)! 



(23) 



This approximation is generally good as long as both Ni 
and N — Ni are large. In our work, N 3> 1, and Ni is 
assumed to be ^ (equal partition) or N Ni 1, so 



A. Equal Partition 

Here for simplicity we set N even (note that the total 
system size is 2N), thus iVi = N 2 = N/2, and \ = 
\{N + y/Nir), then 



X 2 

A Simi 



N + VNn 
and the entanglement entropy becomes, 

E = In | — + -VnN ) 

~lniV-ln4~lnA- 1.38629. 



(27) 



(28) 



Compared with our numerical result, as shown in Fig. [31 
we see that they agree very well, not only for the prefactor 
of the IniV dependence, but the intersection coefficient 
as well. 



B. Unequal Partition 

In this case, 1 <C Ni <C N, we can expand the entropy 
as follows, 



£~ln 



Ni 1 + - 

N 2 



Ni , l>(l-f) 



Ni 



In Ni + In 1 + - 

N 2 



Ni 

InNi-jjr + O 



Nl , l./7T(l-f ) ) (29) 



Ni 
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Numerical Calculation 

— Asymptotic Behavior 
E = lnN-ln4 

— Fit line for Numerical Calculation 
E = 0.99983 lnN-1.0227 
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FIG. 3: (Color online) Comparison of asymptotic behavior 
(black solid line) and exact numerical calculation (scattered 
points) of entanglement entropy E. Exact calculation of E for 
the equal partition case is compared with asymptotic behavior 
we derived in the text up to N — 10000. Results obtained 
from these two methods agree well. Our approximation of the 
asymptotic behavior gives E ~ ln(JV) — 1.38629, compared 
with the best fit line (red dashed line) E = 0.99983 In N - 
1.0227. 




FIG. 4: (Color online) Comparison of asymptotic behavior 
(solid line) and exact numerical calculation (scattered points) 
of entanglement entropy E. As total system size N increases, 
E tends to a constant for a fixed subsystem size Ni, as ex- 
pected. Difference between them increases as the subsystem 
size increases, as the condition 1 <C Ni <C N/2 gets less sat- 
isfied. 



From this expression we see that when the assumed con- 
dition is satisfied the entropy indeed depends only on the 
subsystem size to leading order. To check to what extent 
our approximation is still valid, we plot the expression as 
a function of total system size N with different subsystem 
size Ni , as shown in Fig. 2] This figure shows that when 
TV is large enough the entropy becomes independent of 
N, and the first order term, IniVi dominates except for 
an apparently constant shift from the exact curve. We 
have also plotted the asymptotic behavior as a function 
of Ni while keeping TV fixed, shown in FigO 



FERROMAGNETIC MODEL AND ITS 
ENTANGLEMENT ENTROPY 



In this section we consider a ferromagnetic (FM) spin- 
1/2 model on an arbitrary lattice with N sites, 



H — S JijSi ■ Sj, 



(30) 



with Jij > 0. The ground state is the fully magnetized 
state \SM) with S = N/2 and M = —S, —S + 1, • • • , S, 
and is clearly long-range ordered: (Si ■ Sj) = 1/4. How- 
ever there is a crucial difference between the FM ground 
state and the AFM ground state studied earlier: the FM 
ground state has a finite degeneracy, and thus the sys- 
tem exhibits a non-zero entropy even at zero tempera- 
ture, Eq — log (2 5 + 1) — log(TV + 1), resulting from the 




10000 



FIG. 5: (Color online) Asymptotic behavior of entanglement 
entropy E as function of ATj with N fixed (solid blue line). 
Here N — 10000. Compare it with the scattered points, which 
are an exact numerical calculation of E. We find the asymp- 
totic behavior is accurate until N\ becomes comparable to N. 
We also plot the asymptotic behavior with first order correc- 
tion(dashed light green line) which exhibits the correct tail 
effect when JVi becomes comparable with N. 
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density matrix of the entire system, 

s 

i 

P = 



N 



±— \3M)(SM\ 

M=-S 

^— J2 \SM){SM\ 



N + l 



M=-S 

(\ 
1 

: 1 



(31) 



V 



o 1/ 



M=-S 



In this case the entanglement entropy between two sub- 
systems (1 and 2) is defined in the following manner. We 
first obtain reduced density matrices for subsystems 1 
and 2 by tracing out degrees of freedom in 2 and 1 from 
P- 



1 S 

Pi = tr [2)P = -^—^ J2 tr (2) (\SM){SM\) 

N + 1 ^ 

P^ = tr WP=TFTT E tr (1) (\SM)(SM\) 

A' 

1 5 



PlM, 



(32) 



M=-S 



N- 



M=-S 



and calculate from them the entropy of the subsystems, 
Ei and E 2 . The entanglement entropy is defined as[2fjj 



E=(E 1 +E 2 - Eo)/2. 



(33) 



For the present case E\ and E 2 can be easily obtained 
from the following observations, (i) Because the total 
spin is fully magnetized, so are those in the subsystems: 
Si = iVi/2 and S 2 = N 2 /2. Thus this is a two-spin entan- 
glement problem, (ii) Because the total density matrix p 
is proportional to the identity matrix in the ground state 
subspace, it is invariant under an arbitrary rotation in 
this subspace. (iii) As a result the reduced density ma- 
trix pi is also invariant under rotation in the subspace of 
subsystem 1 with S\ = Ni/2, and is proportional to the 
identity matrix in this subspace. Thus 



Pi = 



N- 



7TT £ 'f*><f«l 

x 1 JV, 



/10 ... 

1 
1 



JVi + 1 



(34) 



and Ei = ln(iVi + 1) (in agreement with Ref. [lj]). Sim- 
ilarly E 2 = ln(N 2 + 1). Thus 

E = [ln(iVi + 1) + ln(7V 2 + 1) - ln(iV + l)]/2. (35) 

We find in both the equal partition (A^i = N 2 = N/2, 
here again we set N even for simplicity) and unequal 
partition (N\ <C N 2 = N — N\) limits, the entropy grows 
logarithmically with subsystem size N\, 



lim E -> -ln(iVi). 



VI. CONCLUDING REMARKS 



(36) 



1/ 



In this paper we have studied antiferromagnetic 
(AFM) and ferromagnetic (FM) spin models with per- 
fect long-range magnetic order in their ground states, 
and calculated the ground state block entanglement en- 
tropy when the system is divided into two subsystems 
(or blocks). In both cases we find the entropy grows log- 
arithmically with block size. In the following we discuss 
the significance of our results. 

First of all, there is entanglement, despite the perfect 
long-range order in the ground state. This is somewhat 
surprising as one might think that in such Hamiltonians 
one can obtain certain properties of the system exactly 
using mean-field theory, in which the ground state is ap- 
proximated by a product state with no entanglement. 
Our results indicate that mean-field approximation is not 
appropriate for entanglement calculation, even if it is "ex- 
act" for other purposes. This point is particularly strik- 
ing for the AFM model, in which the ground state is 
unique. The source of the discrepancy is there still is 
quantum fluctuations even in such a model with super 
long-range interaction, which render the ground state a 
singlet, even though quantum fluctuation does not reduce 
the size of the order parameter. The entanglement is due 
to the quantum fluctuation of the direction of the order 
parameter, which is a collective mode with zero wave- 
vector (or a zero mode); this is missed by any mean- field 
approximation. 

Secondly, the entropy does not obey the area law. The 
reasons for that are different for the two cases we stud- 
ied. In the AFM model, the interaction does not depend 
on distance in the Hamiltonian, thus there is no notion 
of distance or area in this model. In the FM model, 
on the other hand, the ground state is independent of the 
Hamiltonian, as long as all interactions are FM. The fully 
magnetized ground states are invariant under permuta- 
tion of spins, thus there is no notion of distance or area in 
the ground states, even though the terms in Hamiltonian 
can have distance dependence. 

Thirdly, the absence of an area law is special to the 
cases we studied, again each in their own ways. For the 
FM model, it is specific to zero temperature. As soon 
as a finite temperature is turned on, one expects the 
entanglement entropy (or mutual information) to grow 



7 



with the area separating the two subsystems or blocks, 
as long as the interaction is not long-ranged [22| . For 
the AFM model at zero temperature, we do expect an 
area-law contribution to the entropy for short-range or 
even certain power-law long-ranged interaction due to 
quantum fluctuations. This is most easily seen within 
spin-wave approximation, which is a version of mean- 
field theory. Within the spin wave approximation spins 
are mapped onto bosons, and the Hamiltonian is mapped 
onto coupled harmonic oscillators. Detailed recent stud- 
ies have established the area law of entanglement for such 
systems In this regard the super long-range AFM 
model we study here is very special, in that all spin- wave 
degrees of freedom at finite wave- vector disappear, and 
the only degrees of freedom contributing to the ground 
state are zero wave-vector modes represented by Sa and 
S_b [HI ; as a result there is no area-law contribution from 
quantum fluctuations of spin- waves. 

We conclude by speculating that in the more general 
cases that do have an area-law contribution to the en- 
tanglement entropy, as long as long-range spin order is 
present, the logarithmic contribution due to fluctuations 
of the order parameter zero-modes we find here will re- 
main and show up as a sub-leading (yet singular) correc- 
tion to the area law. If that is the case, then conven- 
tional long-range order contributes to the entanglement 
entropy in a way similar to the much subtler topological 
order [3, Q or quantum criticalityQ. 
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APPENDIX A 



6-j Symbol & Racah coefficient: 



9-j Symbol & LS-jj Coupling coefficient: 
'a b c\ 

d e / =((2c+l)(2/+l)(2 fl + l))-* 
K g h ij 

a b c 

x (2h+l)~? d e f ( A2 ) 
_g h i_ 

= ((2c + 1)(2/ + l)(2g + l)(2h + 
x ((ab)c,(de)f;i\(ad)g,(bc)h;i). 



APPENDIX B 

In this Appendix we give a slightly simplified deriva- 
tion of the entanglement entropy for the singlet ground 
state of the infinite-range AFM model defined in Sec. II. 
While this derivation is not as general as that given in the 
main text (which can, in principle, be applied to states 
with nonzero total spin), it has the advantage of clar- 
ifying the reason for the appearance of the Wigner 6j- 
symbol in the final result for the entanglement entropy. 

Following the notation of Sec. II and III, we consider 
the case of 2N spin- 1/2 particles divided equally into two 
sublattices A and B. The spins on each sublattice are 
taken to be fully polarized, so that Sa = N/2 and Sb = 
N/2. In what follows we will use a notation in which, for 
example, the state for which the spins 5,4 and Sb form a 
singlet is represented as (Sa, Sb)o- I n this notation, pairs 
of spins contained within a set of parenthesis form a state 
whose total spin is equal to the subscript labeling the 
parenthesis. Thus, (Sa, Sb)q is a singlet state equivalent 
to the state defined in Eq. @ in Sec. II. (Needless to say, 
for this state to exist it is necessary to have Sa = Sb-) 

Now we consider what happens if these spins arc di- 
vided into two different subsystems labeled 1 and 2. Fol- 
lowing Sec. II, if SAi is the total spin quantum number 
of the A sublattice spins in subsystem i = f,2, and Ssi 
is the total spin quantum number of the B sublattice 
spins in subsystem i ~ 1,2, then the state of the spins 
on the A sublattice can be written (Sa x , Sa 2 )s a and the 
state of the spins on the B sublattice can be written 
(Sbd Sb 2 )s b - The total singlet state for the entire sys- 
tem, \tp), can then be expressed as 



ji h J12 \ _ r _y 1+ j 2+ j 3+ j 
33 J J23' 1 ) 



W(jij 2 Jj 3 ; J12J23) 



= (-i) ji+h+js+J [(2j 12 + i)(2j 23 + 1)]-* 

x (ji, J23; J\Ji2,j3; J)- 



(Al) 



where J12 and J23 refer to the coupling of j\ and j'2 or 
J2 and J3 respectively. 



U') = {(Sa 1 ,Sa 2 ) Sa ,(Sb 1 ,Sb 2 ) Sb ) q - 



(Bl) 



Note that we have not included the m quantum num- 
bers associated with the z-componcnts of Sa and Sb in 
writing the above expression for \tp). This is not neces- 
sary because the requirement that the total spin of the 
two sublattices combine to form a singlet uniquely deter- 
mines the state. It should, of course, always be under- 
stood that when we write Eq. (|B1[) what is really meant 
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(in obvious notation) is, 



7 =L= £ (-ir\(S Al ,S A2 ) SA y, 



(<Sb 15 5b 



(B2) 



Equation (|B2|) effectively (up to irrelevant phase fac- 
tors) gives the Schmidt decomposition (see, e.g., [l|) of 
the state into the two subsystems consisting of the 
A and B sublatticcs. Given such a Schmidt decomposi- 
tion it is straightforward to determine the entanglement 
between these two subsystems. However, here we are in- 
terested not in the trivial entanglement between subsys- 
tems A and B, but the entanglement between subsystems 
1 and 2. To find this we need the Schmidt decomposition 
of into subsystems 1 and 2. 

Note that because is a total spin singlet it follows 
that 



((Sa 1 ,Sa 2 ) Sa , {Sb 1 ,Sb 2 ) Sb ) 
({(S Ai ,S A2 ) Sa ,Sb 2 ) s > s Bi) n , (B3) 



where, for convenience, we have rearranged the order of 
the spins in the second equality. The key point here is 
that the three spins S Al , Sa 2 an d Sb 2 must be in a state 
which is a total spin eigenstate with total spin Sbx ■ If 
this were not the case it would be impossible to form a 
total spin singlet with the remaining spin S b x ■ 

Next we use the Wigner 6j-symbol to express the three 
spin state ((S Al , Sa 2 )s a > Sb 2 )s b as a superposition of 
states of the form (S Al (Sa 2 , Sb 2 ) $ 2 ) s ■ The result is 



((Sa 1 ,Sa 2 ) Sa ,Sb 2 ) Sbi 



(-l) N ^ 7 s 2 (^,(^S B2 ) S2 ) S£ 



where the coefficients are given by 



7*=V(2SA + l)(2S a + l)(f£ |) 



(B4) 



(B5) 



(see Appendix A). 

Rearranging the spins again, and using the fact that 
the total spin of all four spins is 0, we can express the 
resulting four spin basis states as follows, 



iSUi {Sa 2 , Sb 2 )s 2 ) Sbi >&Bi 

= ((Sa 1 ,Sb 1 ) Si ASa 2 ,Sb 2 )s 2 ) (B6) 



where, obviously, Si = 1S2. We therefore conclude that 



= (-!)*£ 7* ((S Ai ,S Bi ) S2 ,(Sa 2 ,Sb 2 ) S2 ) (B7) 



Finally, after writing the m sum explicitly (as in 
Eq. (|B2j) ). we have the desired Schmidt decomposition 
(again, up to irrelevant phases) of the state \ip) into sub- 
systems 1 and 2, 



S m=-S 

<8)\{Sb 2 ,Sa 2 ) s ,-m) 



(B8) 



where 



A Sm = V^Tl(^ %)■ (B9) 



When the values of S A = N/2, S Al = S Bl = N x /2, 
Sa 2 = Sb 2 — (N — N\)/2 are substituted into this ex- 
pression it can be seen to be equivalent to that given in 
Eq. (fT6|) in Sec. III. The derivation of the entanglement 
entropy given in the main text follows. 
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